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i-H Abstract 

o 

Perfect phylogeny consisting of determining the compatibihty of a set of characters is known 
^ to be NP-complete [5l[30]. We propose in this article a conjecture on the necessary and sufficient 

^ conditions of compatibility: Given a set C of r-states full characters, there exists a function /(r) 

such that C is compatible iff every set of /(r) characters of C is compatible. According to 
^ Ii[l0l[i|27l[l2l|2g, /(2) = 2, /(3) = 3 and /(r) > r-1. [25^ conjectured that /(r) = r for any 

VO r > 2. In this paper, we present an example showing that /(4) > 5 and then a closure operation 

for chordal sandwich graphs. The later problem is a common approach of perfect phylogeny. 
This operation can be the ffist step to simplify the problem before solving some particular cases 
/(4), . . . , and determining the function /(r). 

• 

^ Keywords: perfect phylogeny, multi-states characters, chordal completion, triangulation 

^ 1 Introduction 

> 

^ Given an input biological data of a currently-living species set, phylogenetics aims to reconstruct 

y—^ evolutionary history of their ancestors. The evolutionary model of perfect phylogeny is phyloge- 

netic tree, and the data are characters of species. Characters can be morphological, biochemical, 
physiological, behavioural, embryological, or genetic. Each character has several states. Here are 
^ some examples. The character have wings has two states: with wings and without wings. The 

i-H character number of legs has many states: one leg, two legs, four legs, ... These are morphological 

^ characters. For an example of genetic characters, given a set of DNA sequences having a same 

length, if we consider each position on the sequences to be a character, then each character has 4 
states corresponding to 4 bases of DNA as A, T, C, G. 
^ Let £ be a species set, and let c be a character on C. Then, c can be represented by a partition 

of a non-empty subset of C such that each part consists of all species having the same state of c. 
So a set of characters is a set of partitions. A character is said to be trivial if the partition has at 
most one part having more than 1 element. Otherwise, it is non-trivial. If = C, then c is a full 
character, otherwise it is a partial character. If c has at most r parts, then c is a r-states character. 
A binary character is a 2-states full character. 

Example 1 Let C = {a^b^c^d^e^ f^g} and C = {ci, C2, C3, C4} be a collection of characters on C 
such that: 
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ci — ah\cdefj\ghi, C2 — def\ahcghij , C3 — gh\defi, C4 — abcd\ghi. 

It means that ci has 3 states 0, 1, 2 such that the species a, b are in state 0, the species c, d, e, /, j 
are in state 1, and the species g,h,i are in state 2. Similarly for 02,03,04. 

A phylogenetic tree on a species set £ is a tree where each leaf is labehed distinctly by a species 
of/:. 

Definition 1 ^ Let be a r-states character and let T be a phylogenetic tree on C. For i = 
0, . . . , r — 1; denote by Ti(o) the minimal subtree ofT on the leaf set consisting of the species having 
the state i of 0. So, is said to be convex on T iff the subtrees Ti(o) are pairwise vert ex- disjoint. 

A set of characters is compatible iff there exists a phylogenetic tree on which every character is 
convex. So the problem of perfect phylogeny is also commonly known as the character compatibility 
problem. 

Example 2 The set of characters C in Example [I] are compatible because there is a phylogenetic 
tree T in Figure^on which every character is convex. For example, o\ is convex on this tree because 
the subtrees ofT on {a,b}, {o,d,e,f,j} and {g,h,i} are pairwise vertex-disjoint. Similarly, 02, os, 
and 04 are also convex of this tree. 




Figure 1: A phylogenetic tree on which all characters in Example |2] are convex 



Determining the compatibility of a set of characters is NP-complete [5l|30]. In this article, we 
are interested in the necessary and sufficient conditions of compatibility of a set of r-states full 
characters. We propose the following conjecture. 

Conjecture: For any set C of n r-states full characters, there exists a function /(r), which 
does not depend on n, such that C is compatible iff every set of /(r) characters of C is compatible. 

This conjecture is based on the following previous results. According to [TO', |27l [12], /(2) = 2 
and according to [25j, /(3) = 3. Fitch-Meacham examples [8l [9l [271 [25] showed that /(r) > r for 
any r > 2. There are polynomial algorithms of checking compatibility of 3- states full characters 
and 4-states full characters [23j. By [26], if the number of states is restricted, then checking the 
compatibility of C is polynomial. In general, [U (241 [2] showed that there is a polynomial algorithm 
in the number of characters and species, but exponential in the number of states. 
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Some related work 

Existence of perfect phylogeny: Given a set of characters on is there any phylogenetic tree on 
which ah characters are convex? It is easy to determine whether a cohection of binary characters is 
compatible. However the problem is NP-complete even for 2-states characters [5;i30j. There exists 
effective, practical approaches for 2-states characters [19], and a new approach for this problem is 
proposed in [14j. 

Quartet problem: A minimal non-trivial character is a 2-states character such that each state 
contains exactly two species. Such character is called a quartet. As stated in the previous paragraph, 
the problem of compatibility of a set of quartets is NP-complete [5^.30j. However, there are some 
particular cases that the problem is polynomial [3j, see [28j for details. 

Define a tree by characters: a set of characters define a tree iff there is not any other tree on 
which these characters are convex. showed that a set of 3 characters are not sufficient to define 
a tree but a set of 5 characters are. Later, [22j showed that for any tree, there exist at most 4 
characters which define this tree. Hence, 4 is the optimal value. For the problem of whether a set of 
characters defines a tree, this is recently proved to be NP-hard [20]. The author solved the quartet 
challenge proposed by [28j as follows: given a phylogenetic tree on C, and a quartet set Q which is 
convex on this tree. Is there any other tree on which this quartet set is also convex? Despite the 
NP-hardness, there is a polynomial algorithm for this problem when — |Q| = 3 

Maximum parsimony: When there is no perfect phylogeny that can be inferred from data, it is 
desirable to find a model that minimize the number of reverse and convergent transitions. That is 
the problem of maximum parsimony. 

Perfect phylogeny with recombination: When the characters set are not tree-representable, it 
is also interesting to construct a model that can represent phylogenies. The model used here is 
recombination networks. Introduced by [21j, then intensive work have been done on this problem, 

including i3ii [ig [ig [13 [13 [E] . 
2 Preliminaries 

A very popular approach of perfect phylogeny is using chordal completion of vertex-coloured graphs, 
or equivalently chordal sandwich graph. 

Definition 2 Let C = {xi, . . . , x^} be a species set and let C = {ci, . . . , c^} be a set of characters 
on C. Each Ci is a partition of a subset of C. The partition intersection graph G — {V^E) of 
C is constructed as follows: 

- Each character of C is associated with a different colour. 

- Each vertex of V corresponds to a state of a character of C. This vertex is then coloured by 
the colour of the character. 

- There is an edge between 2 vertices if the 2 corresponding states of the 2 characters have at 
least a common species. 

In our figures, in stead of colouring the vertices, we include the name of the characters in the 
labels of the vertices. 

Example 3 Let consider the character set in Example^- C = {ci, C2, C3, C4} where 
ci = ab\cdefj\ghi = ci,o|ci,i|ci,2, 
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C2 = def\abcghij = C2,o|c2,i, 
C3 = gh\defi = C3,o|c3,i; 
C4 = abcd\ghi = C4q|c4^i. 

E'ac/i vertex Cij represents the state j of character i. The partition intersection graph of C is in 
Figure 

^1,1 ^,1 ^1,2 



C2,o C3.1 C4-, 



Figure 2: An example of partition intersection graph 



A graph G is chordal if every cycle of length > 4 contains at least a chord. A chordal completion 
of G is a chordal graph G' — {V^E') such that E C E\ This completion is minimal iff when we 
remove any edge in E' \ E^ the resulting graph is not chordal. 

Given a vertex-coloured graph G, a proper chordal completion of G is a chordal graph 
G' — (y, E') such that E <Z E' and E' does not contain any edge connecting two vertices of the 
same colour. 

Theorem 1 |^ \3U^ A set of characters C is compatible iff its partition intersection graph has 
at least a proper chordal completion. 

The set of characters in Example |3] is compatible as showed in Example [2j Its partition inter- 
section graph (Figure [2]) has indeed a proper choral completion which is itself. 

Proper chordal completion of vertex-coloured graph can be stated equivalently under the form 
of sandwich problems. 

Definition 3 Given a graph G = (V, E, F) where F is a set of pairs of vertices of G such that 

EnF^d). 

If there is a graph Gs = {V, Es) such that E C Es C E x E \ F and Gs satisfies property li, 
then Gs is called a U-sandwich graph of G. 

See [llj for some problems and results on graph sandwich problems. 

It is easy to see that a chordal completion of a vertex-coloured graph G = {V, E) is proper iff it 
is a chordal-sandwich graph of G = (V, F) where F is the set of pairs of vertices having a same 
colour. So, by considering the set F, we can ignore the colours of the initial graph. We also call 
a chordal-sandwich graph of G a proper chordal completion of G, i.e. a chordal completion of G 
without using any pair of vertices in F. 

3 Our contributions 

Fitch-Meacham examples were first introduced in [HI [9] , then later generalized in [27] and formally 
proved in [25], showed that /(r) > r for any r > 2. [25j conjectured that for any r, there is a 
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perfect phylogeny on r-state characters if and only if there is one for every subset of r characters, 
i.e. /(r) = r for any r > 2. However, we have an example in Section [4] showed that /(4) > 5. It 
improves the lower bound of Fitch-Meacham examples and shows that the conjecture in [25j is not 
true. After that, in Section [5) we propose a closure chordal sandwich graph operation such that 
the obtained graph has a stronger structure. As a consequence, one can suppose that the input 
graph has such a structure, which can facilitate settling the conjecture. 

4 An example of 4-States Characters 

We present here an example of a set of 4-states characters which is not compatible, but every 4 
characters of this set are compatible. 




Figure 4: The induced subgraph of G on 4 colours 
Figure 3: Graph G a, c, d, e. This graph is chordal. 

Let C be the following set of characters: 
a = {x,u}\{z,t}\{y}\{v} = ao|ai|a2|a3 
6={x,y}|{t,i;}|{z}|M = 6o|6i|62|63 
c = {y,z}\{u,v}\{x}\{t} = co|ci|c2|c3 
d = {x,u}\{y,z}\{t}\{v} = d^\di\d2\d^ 
e = = eo|ei|e2|e3 

Each character has 4 states that we denote for example by ao, ai, a2, for the character a. The 
partition intersection graph G associated to C is in Figure [3} 




Figure 5: The induced subgraph of G on 4 Figure 6: The induced subgraph of G on 4 
colours a, 6, c, d and its chordal completion colours 6, c, d, e and its chordal completion 
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G does not accept any proper chordal completion. Indeed, if we consider only the induced 
subgraph of G on 4 colours a, 6, c, d and triangulate it, then there is a unique way to do that by 
connecting (ai, 60)5 (^o, ci) and (ci,ai) (Figure [5]). Similarly, if we consider the induced subgraph 
of G on 4 colours a, 6, c, e, there is also a unique way to triangulate it by connecting (ao, cq) 
and (co,ao). So, to triangulate G, the cycle (ao^iai^o, ag) is forced to be created. However, this 
cycle does not have any proper chordal completion. In other words, G does not have any proper 
chordal completion. So, C is not compatible. 

However, as we see in Figure |4j the induced subgraph of G on 4 colours a, c, d, e is chordal. 
The induced subgraph of G on 4 colours a, 6, c, d has a proper chordal completion (Figure [5]) and 
similarly for the induced subgraph of G on 4 colours a,b,c,e due to the symmetry. The induced 
subgraph of G on 4 colours 6, c, d, e also has a proper chordal completion (Figure |6]) and similarly 
for the induced subgraph of G on 4 colours a, 6, d, e. 

It means that every 4 characters of C are compatible but the whole set C is not compatible. 

5 A closure chordal sandwich graph operation 

Given a graph G = (V, E, F) where E H F = 9. Let u, v be two vertices of G such that {u, v) E. 

{u^ v) is a forbidden edge if it is not included in any minimal proper chordal completion of G. So 
(u, v) is forbidden if either {u, v) ^ F oi if by connecting them, the resulting graph does not have 
any proper chordal completion. The forbidden edges are presented by dotted lines in our figures. 

(u, v) is a forced edge if it is contained in every proper chordal completion of G. So if there is 
a cycle in G which has a unique proper chordal completion, then the edges used to complete this 
cycle are forced. 

Note that by adding any forced edge into E or any forbidden edge into F, we do not lose any 
proper chordal completion. 

A cycle G of G is forbidden if every chordal completion of G contains at least an edge in F. 
So if G has a proper chordal completion then it does not have any forbidden cycle. The converse 
is not always true. For example, see the graph G = (V^E^F) in Figure [s] where F consists of the 
pairs of vertices having a same colour. This graph has 3 chordless cycles and each one can be 
chordally completed without using any edges in F. However, G does not admit any proper chordal 
completion. 

Example 4 In Figure we have a cycle of size 5 on 3 colours a, 6, c. The set F consists of 
(ao,ai) and (bo, hi). One can deduce that (ai,6i) is forbidden since by connecting them we have 
the forbidden cycle (aoboaibi). We deduce furthermore that (co,ao) and (co,6o) ^^^e forced because 
the unique way to properly chordally complete this cycle is connecting them. 




Figure 7: (ao, ai), (60, ^1), (^i, 61) are forbidden edges. Figure 8: A graph which does not have 
(co,ao), (co,6o) are forced edges. any proper chordal completion 
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Consider the graph in Figure where F is the set of pairs of vertices having a same colour. 
Similarly to the previous example on the cycle aobicoaibo, we deduce that coao and cq^o forced. 
So, the cycle 6oCo^2Ci is forced to be presented in any proper chordal completion of this graph. 
However, this cycle is forbidden. Therefore, this graph does not have any proper chordal completion. 

Observation 1 Given a cycle C = (ui . . .Uk^ui), then: 

(i) for any i G every chordal completion of C must contain either i/^+i) or 

[ui^ Uj) for a certain j different from z, z — 1, z + 1. 

(a) every chordal completion of C must contain a chord (ui-i^Ui^i) for a certain i. 

Lemma 1 (Detecting forbidden edges and forced edges) Let G — {V^ E, F) be a graph where 
E n F = 9 and {u, v) be two vertices of G: 

1) If there is a chordless path (uti . . . t^v) such that for any i — 1^ . . .k, either {u^ ti) or {v^ ti) 
is forbidden, then {u, v) is also forbidden. 

2) Suppose that (u^v) G E. If there is a chordless cycle c = (uwti . . .tkV^u) such that for any 
i — 1^ . . .k, either (u, ti) or {v, ti) is forbidden, then {v, w) is a forced edge. 




uv uv uv uv 

(a) {u,v) is for- (b) {v,w) is a (a) a /(i^, path, (iz, is for- (b) a g{u,v,w) cycle, 

bidden forced edge bidden. {w, v) is forced. 



Figure 9: Lemma [T] Figure 10: Corollary [T] 

Proof: 1) By connecting (u^v), we obtain the chordless cycle C = {uti . . .t^v^u). We will prove 
the cycle C is forbidden. By Observation [l] (i), to chordally complete C, we must connect either 
{ti^v) or {u^ti) for a certain i — 2, . . . , fc. However, {ti^v) is forbidden because (u^ti) is an edge 
of G and by the assumption, either {u,ti) or {v,ti) must be forbidden. So, we must connect an 
edge {u^ti), which must not be a forbidden edge. We deduce that {v^ti) is forbidden. The created 
subcycle {uuiUi^i . . . UkV ^ u) has the same property as C. So, by using the same argument, to 
chordally complete this cycle, we must connect an edge {u, uj) where i < j < k. So, the size of the 
considering cycle decreases strictly each time applying this argument. Then there will be a moment 
that we obtain a cycle which is forbidden, in other words the cycle C is forbidden. 

2) We will prove that, every proper chordal completion of C must contain the chord {v^w). 
Suppose that there is a proper chordal completion of this cycle which does not contain {v,w). By 
Observation [l] (i), this completion must contain {u, ti) for a certain z = 1, . . . , A:. We obtain then a 
subcycle (utiti^i . . . tkV, u). The path (utiti^i . . . tkv) satisfies the condition of Claim 1, so (u, v) is 
a forbidden edge. However, (u^ v) is an edge of G. It means that this subcycle does not have any 
proper chordal completion. In other words, C does not have any proper chordal completion which 
does not contain (u^v). So (v^w) is presented in any chordal completion of G, i.e. it is a forced 
edge. □ 
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Corollary 1 Let a graph G — {V^ E, F), denote by F{u) the set of vertices u' such that {u, u') G F . 
Then, for any pair of vertices {u^ v): 

1) If {u,v) is not an edge of G and there is a chordless path {uti . . .tkv) such that for any 
z = 1, . . . , fc; G F(u) U F{v) then (u^ v) is forbidden. We call such a path a f(u^ v) path. 

2) If (u^v) is an edge of G and there is a chordless cycle (uwti . . .t^v^u) such that w F(v) 
and for any z = 1, . . . , fc^ G U Fy, then {v, w) is a forced edge. We call such a cycle a g{u, v, w) 
cycle. 

Denote by N{u) the set of neighbour vertices of u. 



Data: A graph G = {V, E, F) 
Result: Closure{G) 

1 For any u calculate N{u) and F{u) — {v\ {u,v) G F}; 

2 flag = true; 

3 while (flag) do 



4 
5 
6 
7 
8 
9 

10 
11 

12 
13 
14 



flag = false; 

for (any pair of vertices (u^ v) ) do 
if (there is a f{u^v) path) then 
if (u N{v) ) then 

Add u to F{v)^ and v to F{u)-^ 
flag=true; 

else 

^ G does not have any proper chordal completion. Exit] 

if (u G N{v)) A (there is a g{u,v,w) cycle) then 
Add V to N{w)^ and w to N{v); 
flag=true; 



15 return = {V,E\F') where E' = {(^,^^)| u G N{v)} and F' = {(^,^^)| u G F{v)}. 



Algorithm 1: A closure chordal-sandwich graph operation 



Theorem 2 Algorithm^ takes time 0(n^(n + m)). Let G' = (V^E'^F') = Closure{G), then any 
proper chordal completion of G' is a proper chordal completion of G and vice-versa. Moreover, G' 
satisfies the following properties: 

1. For any {u,v) ^ E' {J F' , if we connect (u,v) then for any created chordless cycle C which 
has (u^v) as an edge, C has at least a chordal completion without using any pair of vertices in F' . 

2. Any chordless cycle of G' has at least two chordal completions without using any pair of 
vertices in F^ 



Proof: Correctness: 

According to Corollary [l] the loop for recognizes all pairs of vertices [u^ v) which satisfy the 
conditions in the corollary to detect the forbidden edges and forced edges. So the pairs of vertices 
added in F' at line 8 are the forbidden edges, i.e. the edges which are not included in any proper 
chordal completion of G. And the pairs of vertices added in E^ at line 13 are forced edges, i.e. the 
edges which are included in every proper chordal completion of G. Therefore, any proper chordal 
completion of the obtained graph is a proper chordal completion of G and vice versa. 
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Moreover, if {u,v) is an edge and there exists a f{u,v) path then according to the proof of 
Lemma [!} the chordless cycle consisting of f{u^v) and {u^v) is forbidden. Hence, G does not have 
any proper chordal completion (line 11). This process stops when there is no more forbidden edges 
or forced edges detected. So, in G\ for any pair of vertices (u^v), there is no f(u^v) path; and if 
(u^v) G E\ then there is no g(u^v^w) cycle. We will prove that satisfies the two properties of 
the theorem. 



1\ 












wl / 


v. 



Figure 11: Proof of Theorem [2] 



1) The first property: we prove by induction on the size of cycles. 

Let (u^v) be a pair of vertices not in U F^ By connecting {u^v), let C be a created cycle 
which contains (u^v). We will prove that C admits at least a chordal completion without using 
any edge in F^ 

For the case \C\ = 4, let C = (uxyv^u). By the assumption, G^ does not contain any f{u^v) 
path, i.e. (uxyv) is not a f(u,v) path. It means that either {v,x) or (u^y) F^ So, we can 
complete C by connecting either (v^x) or (u^y). 

Suppose that C has at least a chordal completion if |C| < A:. 

For the case |C| = + 1, let C = {uti . . . t^-iv, u), so there is not any f{u^ v) path in G^ Then, 
there exists at least an i G {1, . . . , A: — 1} such that (u, U), {v, U) F^ By the induction hypothesis, 
if we connect (u^ti), i^^^ti), the two subcycles (uti . . .ti^u) and {vti . . .tk^v) have proper chordal 
completions without using any pair of vertices in F' because both these two cycles have size smaller 



than + 1 (Figure 11(a)). Completing these two subcycles gives a proper chordal completion for 
C. So, C admits at least one chordal completion, the one which contains {u^ti) and (v^U). 

2) The second property: Suppose that there is a chordless cycle C which admits a chordal 
completion without connecting any pair of vertices in F^ By Observation [l] (ii), this chordal 
completion must contain at least a triangle (u^v^w) such that {u,v), {u,w) are edges of C. Let 
C = {uwti . . .tkV^u), so C is not a g{u^v^w) cycle because otherwise (v^w) is a forced edge and 
it must have been connected by the algorithm, a contradiction with the fact that C is chordless. 



So, there is a U such that {u^U), {v^U) (Figure 1 (b) ). Using the first property, if we connect 



(u^ti) and (v^ti)^ then we obtain two subcycles which have proper chordal completions without 
connecting any pair of vertices in F^ That implies another chordal completion of C containing 
{u,ti) and {v,ti). This chordal completion does not contain {v,w), so it is different with the initial 
one. In other words, C has at least 2 distinct proper chordal completions without using any pair 
of vertices in F^ 

Complexity: 

- Calculating N(u) and F{u) for any vertex u in line 1 is done in times O(n^). 
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- The loop while: For each iteration, there is at least a pair of vertices {u,v) whose nature is 
modified, i.e (ix, v) becomes either a forbidden edge or a forced edge. Once it is modified, it will 
not be modified afterwards. The loop stops when there is no more modification on any pair of 
vertices. So, the number of iterations of this loop is bounded by the number of pairs of vertices, i.e 
by 0(n2). 

- The loop for: there are O(n^) pairs of vertices {u,v). So there are 0{n?) iterations. In each 
iteration: 

• Checking if there is a simple path f{u,v) can be done in linear time: We proceed a dfs 
starting at u such that the visited vertices are in F(u) U F(v) \ N(u). If we meet a vertex in N(v)^ 
then there is a f{u^v) path. Otherwise, there is no such path. 

• Checking if there is a g{u^v^w) cycle can also be done in linear time: We proceed a dfs 
from u such that the first visited vertices is not in N{v) U F{v)^ and the remaining visited vertices 
are in F{u) U F{v) \ N{u). If we meet a vertex in N{v) then we have a g{u^ w) cycle. Otherwise, 
there is no such cycle. 

So, the total complexity is 0{n^{n + m)) where n is the number of vertices of and m is the 
number of edges of the obtained graph. □ 

With this operation, one can deduce for example that the input graph does not contain any 
cycle as in Figure [7| because this cycle must be already triangulate by the operation. Or the input 
graph can not contain an induced subgraph as in Figure [8) because the operation can show that in 
this case the graph does not have any proper chordal completion. 

Corollary 2 Without loss of generality, one can suppose that any cycle of a graph G — {V^ F) 
has at least two proper chordal completions; and by connecting any pair of vertices not in E U F, 
every created cycle has at least one proper chordal completion. 

6 Conclusion 

Our example in Section [4] showed that /(4) > 5. We suggest that /(r) > r + 1 for any r > 4. So, 
a further work is to generalize our example, or to find other examples supporting this suggestion. 
Another problem is proving that /(r) exists by determining an upper bound function of r for /(r). 
It means that we must find a function F{r) such that if every set F(r) characters of C is compatible 
then C is compatible. A harder question is determining /(r) for r > 4. Our closure operation for 
chordal sandwich graphs can help to simplify these problems. 
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